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TECHNICAL NOTE 3609 


LINEARIZED LIFTING-SURFACE AND LIFTING-LINE EVALUATIONS 
OF SIDEWASH BEHIND ROLLING TRIANGULAR 
WINGS AT SUPERSONIC SPEEDS 
By Percy J. Bobbitt 


SUMMARY 


The lifting-surface side-wash behind rolling triangular wings has 
been derived for a range of supersonic Mach numbers for which the wing 
leading edges remain swept behind the Mach cone emanating from the wing 
apex. Variations of the sidewash with longitudinal distance in the 
vertical plane of symmetry are presented in graphical form. 

An approximate expression for the sidewash has been developed by 
means of an approach using a horse shoe -vortex approximate -lifting -line 
theory. By use of this approximate expression, sidewash may be computed 
for wings of arbitrary plan form and span loading. A comparison of the 
sidewash computed by lifting-surface and lifting-line expressions for 
the triangular wing showed good agreement except in the vicinity of the 
trailing edge when the leading edge approached the sonic condition. 

An illustrative calculation has been made of the force induced by 
the wing sidewash on a vertical tail located in various longitudinal 
positions . 


INTRODUCTION 


In order to make reliable estimates of the total forces and moments 
acting on an aircraft, accurate evaluations are required of the loadings 
on the individual isolated components and of the interference effects 
between components. Although considerable effort has been expended in 
recent years to supply much of this needed information for the supersonic 
speed range, many important problems remain. Among these is the induced 
effect of the wing flow field or, more precisely, the wing sidewash on the 
vertical tail. The only specific numerical results of this nature obtained 
to date have been for the angle -of -attack motion. In the vertical plane of 
symmetry for this case, however, the sidewash is zero and tail surfaces 
located in this plane are unaffected. This is not the’ situation for the 
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rolling, yawing, and sideslipping motions where the sidewash in the ver- 
tical plane of s ymme try is finite and the load induced on the vertical 
tail can be appreciable. Evaluation of the sidewash for these motions 
would, therefore, be Imp ortant in the prediction of the lateral stability 
of supersonic aircraft . 

The present paper presents the derivation of the sidewash behind 
steady rolling, triangular wings with subsonic leading edges. Both lifting- 
surface and lift ing -line methods, previously applied primarily to determine 
downwash, are utilized and comparisons are made of the sidewash computed 
by the two methods in order to give an indication of the worth of the more 
easily obtainable lift ing - line results . The lifting-surface sidewash is 
determined by using the doublet-distribution method of reference 1, and 
the lifting-line values are obtained by use of the lifting -line approach 
given in reference 2. 

An illustrative calculation using the derived sidewash is made of the 
force induced on a half -delta tail operating behind a rolling triangular 
wing, and this force is contrasted to the force that would act on the tail 
if it were rolling in the undisturbed stream. 

Thp material presented in this paper was submitted to the University 
of Virginia as a thesis in partial fulfillment of the requirements for 
the degree of Master of Science in Aeronautical Engineering. 


SYMBOLS 


The positive directions of forces, moments, and velocities are shown 
in figure 1. 

x, y, z Cartesian coordinates of field point 

x l, yi> z i Cartesian coordinates of doublet or line-vortex position 


u, v, w 
A 
b 
^n 


perturbation velocities along x-, y-, and z-axls, respectively 

t 

wing aspect ratio, b^/ s 
wing span 


yawing-moment coefficient. 


Yawing moment 
qSb 
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dc r 

5 * 

2V 


p — >0 


C Y 


side -force coefficient , 


Side force 
ciS 




2V 


c 

d 

h 


P — ^0 

wing root chord 

distance from wing trailing edge to a point downstream 
displacement of vortex sheet "below wing trailing edge 


hi, hg limits of y^-integration 

i variable index used in s umma tions and as subscript 


k 


k l = 


constant 


= ^/( X ~ x l ) 2 - P 2z ' 

e o x l 


%° 


__ ^/(x - c) 2 - P 2 z 2 


0 o c 


%o' = c / 1 - k l/ 


**2 “ 


e o x l 


7(x - x i) 2 - p 


2 2 
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M free-stream Mach number, V/Velocity of sound in free stream 

m slope of lifting line 

P static pressure 

AP = P Z - P u 

p angular velocity of roll, radians/sec 

q free-stream dynamic pressure, | PV 2 

S wing area 

Au s = % - 

V free-stream velocity 

vp sidewash Induced by doublets distributed over plan form 

v^ sidewash induced by doublets distributed over wake 

X = x - xp 
x i = x - Xl 



x 

x o = C 

X 1 

x l,° = c 
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y = y - y x 

y i = y - Yi 

Y -II 

i ^° D/2 

y = -4- 
0 b 72 


z 


o " 


b/2 


a angle of attack, radians 

p = \Jyfi - 1 

T circulation at any spanwise station 

6 angle of downwash measured in xz -plane, "between trailing 

vortex sheet- and axis parallel to free-stream direction, 
radians 

0 O = ^ for triangular wingsj 

A angle through which vortex sheet rotates in moving from wing 

trailing edge, to a point which is d distance downstream 

p density of free-stream air 

(/> perturbation velocity potential 

A 0s - 0u - h 

n - 6 0 C - 1 

pz Zq 


E 


complete elliptic integral of second kind with modulus k, 
P L \/l - k 2 s 2 


ds 


j 0 y/l - s 2 


K 


complete elliptic integral of first kind with modulus k, 
fU 

ds 


UO Jl _ k 2 s 2 \|l - s 2 
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E(t,k) incomplete elliptic integral of second kind with argument t 

ds 


and modulus k, 


\J± - k 2 s 2 \jl - s 2 


F(t,k) incomplete elliptic integral of first kind with argument t 

Pt 


and modulus k. 


\/l - k 2 s 2 ds 

Jo - B 2 


f x - Jq q 2 '^ + p 2 z 2 (i-e 0 2 ) 


1 - 9r 


_ x o - 9 o / x o 2 + z o 2 ( 1 ~ 9 o? 


f O = 2 

1 - 0 O 2 


G ( k 2,o^^ - | + (* 2,0 " E 2,o) F ^ ot ^ *2,0^ ~ k 2>0 j 


2-0 


o( 9 o) - °-E'(e 0 ) ^-f’( 9 °) 


1-0 


1 - 9f 


P a a(0 o ) 

L.E. leading edge 

T.E. trailing edge 

sign denoting finite part of integral 
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Subscripts : 

D conditions in region D (fig- 2) 

E conditions in region E (fig. 2) 

l pertaj.ning to lower side of surface 

P plan form 

s conditions on surface of discontinuity (at Z]_ = oj 

u pertal.ni.ng to upper side of surface 

W wake 

The subscripts 1,2, l,o, and 2,o on the elliptic functions E 
and K indicate the modulus of the elliptic function; that is, 

% = K (f>*i) 

E 2,o ' E (f,*2,o) 


ANALYSIS 
General Remarks 


The problem to be considered herein is that of determining the 
perturbation sidewash velocity behind a rolling triangular wing for a 
range of supersonic Mach numbers for which the leading edges of the 
wing are subsonic. The analysis is based on an application of linearized 
supersonic-flow theory and, hence, the results obtained will be valid 
within the limitations of linear theory. 

In the analysis which follows, several assumptions are made con- 
cerning the trailing vortex sheet. These assumptions are that the vortex 
sheet must remain flat behind the wing and that the rotation of the 
vortex sheet is small enough to be neglected. In addition, the nonre- 
strictive stipulation is made that the roll i ng wings be at zero angle of 
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attack. Further discussion of these points will he found in the section 
entitled "Discussion. " 

In flight, a steady rolling motion will usually he maintained hy 
differentially deflected ailerons that create a sidewash opposed to the 
wing sidewash. Calc u l a tion of aileron sidewash, which may he of the 
same order of magnitude as the wing sidewash, will not he considered in 
the present paper. 


Boundary Conditions 

The houndary conditions for the proposed problem may he prescribed 
on the z = 0 plane and are similar to those given for the angle -of -attack 
motion in reference 1. 

The downwash boundary condition on the rolling wing is 


w* = 



In order to analyze the quasi-steady rolling problem hy use of steady- 
flow theory, the rolling wing is considered fixed in approximately the 
z = 0 plane hut twisted linearly in the spanwise direction. Only small 
linear twists are allowable, however, in order not to violate the assump- 
tions of small-perturbation linearized theory; hence, the rate of roll is 
necessarily small (approaching zero). 

Pressures on the wing and pressure differences across the wing surface 
are finite and, for a great variety of plan forms, have already been 
obtained. (See, for example, refs. 3 and k.) Off the wing and in the 
plane, of the wing, the pressure, and hence the pressure difference, must 
he zero. 

In the z = 0 plane, the local pressure difference is directly pro- 
portional to the streamwise component of the perturbation velocity and is 
given simply as 


^(*1^1) _ 2 Au s( x i^i) 

V 


q. 


(1) 
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By consideration of the relationship Between the perturbation velocity 
potential and the streamwise velocity component; that is , 




( 2 ) 


an expression giv ing the jump in velocity potential across the xy-plane 
in terms of the local pressure difference may be written as r 


^s( x i^i) 



^ X 1 L 
l L • 


C 1 APfx^yJ 


dx 


<3o 


(3) 


Since, from equations (l) and (2), 


= 2 ^ ^s( x l^l) 

q V chcp 


w 


and since 0 B ^i>yiJ is an odd function in z, it is clear that, beyond 
the trailing edge, A0 S must be independent of xp to satisfy the zero- 

pressure condition in the wake . The integration indicated in equation ( 3 ) 
should, therefore, be made from the wing leading edge to the trail in g edge 
to obtain A0 a in the wake. 


Solution to Boundary -Value Problem 

The linearized partial-differential equation which the perturbation 
velocity potential must satisfy in supersonic flow is 


2 5^0 c ^0 ^0 _ 

3x^ dy^ dz^ 


0 


(5) 
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For the problem "being considered herein, the solution to equation (5) may 
he written as 



This expression represents the potential in space due to a distribution 
of doublets in the xy-plane with strengths that are governed by the 
potential jump across the z = 0 plane. 

The symbol [ indicates that the finite parts of the infinite 

integrals are to be taken when they appear. 


Application to Triangular Wings 

The loading over a rolling delta wing with subsonic leading edges 
has been found in reference 4 to be 



%x 1 e 0 2 gy 1 

VP 2 g(q 0 ) 


( 7 ) 


where 


2 

1 





G ( e o) - 
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and 

E '( e o) - e (| >fi - e o 2 ) 

r '( e o) - - »?) 

From equations ( 3 ) and ( 7 ) the potential jump across the wing surface is 


A0 s (x lJ y 1 ) = Hp yi JQq 2 *! 2 - P 2 Yi 2 (8) 

* 

and in the wake 

^ s (yJ = Hpy 1 ^e 0 2 c 2 - p 2 y x 2 (9) 

where 


= 

p 2 c(0 o ) 


( 10 ) 


The velocity potential in space may now he written as the sum of the 
two expressions (see eq.. ( 6 )) 


zH{3 2 

rr 

Py l ~ P 2y i 2<3x i 

2jt 

JJ 1 

Plan form 1 

[(x - xp) 2 - p 2 (y - yp ) 2 - p 2 z 2 5//2 


( 11 ) 
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and 


zH3 2 

rr 

pyi \l 0 o 2c2 - <3yp 

2it 

JJ r. 

Wake I f: 

s - Xp ) 2 - p 2 (y - yj 2 - P 2 z 2 5//2 



As previously stated, the primary purpose of this paper is the 
determination of the velocity perturbated behind the wing parallel to 
the y-axis (or the sidewash) . This flow velocity may be obtained by 
taking the partial derivative of the velocity potential with respect 
to y, or ; 


v 


__ cffi 
dy 


With 


0 = 0 P + 0^ 


the sidewash in the xz -plane from equations (ll) and (12) will be given 
by the sum of 





sH3 2 S 


- lim _ 

y — >0 2jt oy 


2 y 2 
X 1 


P 2 y 1 2 dx 1 dy x 


P y l & 

Plan form ^x - x^ 2 - g 2 ^y - 2 - p 2 z.2j "^ 2 


( 13 ) 
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and 



-lim *SLA\ 

y — ^0 2rt dy 



pyi \J e o (l) 2 c 2 - P 2 ^^! foi 

(x - xj 2 - p 2 (y - yi) 2 - p 2 z2j 5 / 2 


(3A-) 


Subsequently it -will be convenient for computational purposes to 
derive expressions for the sidewash' which hap been nondimensionalized by 
pb/2 so that 


v = V P + V W 
pb/2 pb/2 pb/2 


When v is written in a slightly different, though equivalent, form 


v/V 

pb/2V 


it can be recognized that the nondimens ional sidewash parameter may 
be defined as the induced angle of sidewash per unit wing-tip helix 
angle pb/2V. 

The rest of this section is devoted to the evaluation of equa- 
tions (13) and (14) at points in the following two regions of the 
xz-plane (see fig. 2): 

(l) The region lying between the Mach lines emanating from the wing 

trailing edge and the line of intersection of the two cones from the 

trailing -edge tips . 
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(2) The region which extends from the line of intersection of the 
two cones from the trailing-edge tips downstream to infinity. 

These two regions are denoted, as in reference 1, "by E and D, 
respectively. The contributions of the doublets distributed over the 
plan form and the wake to the sidewash in regions E and D are con- 
sidered separately. 

Sidewash due to doublets distributed on plan form in region D . - 
In region D the sidewash contributed, by the doublets distributed on the 
plan form is 


d0P,D = 
dy 


6 0*1 


= - lim 

y — >01 2it 


a rr 

*y Jo 


P pyj s] Q q 2x i 2 - P^! 2 ^ dx 1 

9 o x l ^x - x^ 2 - p 2 ^y - y^ 2 - p 2 z 2 ^/ 2 J 


( 15 ) 


In order to facilitate the integrations involved in determining v p it 

r ,u 

is convenient to carry out the differentiation and limiting processes 
first. This procedure gives 




_ ^hp 2 r c 
* Jo 


9 o x l , 

f p i3^y L £ Je Q g ^i g - s 2 yi 2 flyi 

1 [(* - Xl )2 . pV - 

— _ 


(16) 


When the following substitutions are made. 


2 _ P^l 2 


s = 


fl 2 V 2 
e o X 1 


l ^ 2 = go 2 ^ 2 


(x - X x j 2 - p 2 z 2 
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V P,D "becomes 


= im 2 *2^1 f 1 s 2 y/ L _ s 2^ 

* JO 9 o x l Jo (l - k 2 2 s 2 ) 5 / 2 


(17) 


which, by the use of partial fractions, . may be written in the more 
amenable form 


'P,D 


3zHp <: 




P c k^dxp 

-P - a 

Jo 6 0*1 

_ Jo 



(18) 


The integrals in equation (18) may be reduced to standard elliptic 
forms by use of the Jacobian transformation, s = sn u (refs. 5 and 6), 
and readily integrated to give 



+ 2 


-Eo dXn 


1 - 


(19) 


By replacing H by its equi valent and nondimensionalizing xp, x, and 
z, equation (19) becomes 
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v P,D -2 

pb/2 




( 20 ) 


where 

- : °° 2x i 2 . e o 2x i,° 

(x - X1) 2 - u a z 2 (x Q - *l,o ) 2 - e o £z o 2 

Sidewash due to doublets distributed on plan form In region E .- 

The portion of the wing area over which the integration in equation (13) 
is to be performed is different for each position of the field point in 
region E. This fact is evidenced by the appearance of the field-point 
coordinates x, y, and z in the limits of integration. It is expedient 
in determining Vp -g to follow the same procedure used in determining 

Vp p (differentiating before integrating) . This is allowable since it can 
be shown that the expression for the potential 0p^ E can be differentiated 

with respect to y without regard to the variable limitB when the evalua- 
tion of the derivative is made at y = 0. 

Differentiating </L with respect to y and then setting y equal 

Jr,* 

to zero yields 



( 21 ) 
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where 


f = 


x 


- /Qq 2 * 2 + P 2z2 ( x - e 0 g) 


1 


- e 


2 

o 


The y^-integration in the first double integral of equation (21) is 
identical to the y-]_-integration in equation (l6); hence, only the second 
term of equation (21) remains to be considered. 

In the integration 


let 




/e 0 2x i 2 2 

p _ 

( x - X 1 

l - e 2 * 2 . . si 

5/2 

L e 2 


. yi 2 

0 


where 


2 _ ( x ~ x l) £ - P 2g 


Q 2y 2 
°0 X 1 
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This substituion results in the expression 


P 

6 o x l k l 2 



t 2 \j± - kT 2 t 2 at 


(X - t?) 


5/2 


( 22 ) 


which, except for the use of the finite -pant concept, may he integrated 
in a manner similar to that used for equation (17)- Performing the 
integration in expression (22) allows Vp ^ (eq. (21)) to he written as 



Equation (23) for Vp^-g does not lend itself readily to numerical cal- 
culations because, at the point x^ = f, k-^ and kg become equal to 

1 and give rise to first-order infinities in the integrands. Appendix A 
shows how this difficulty is eliminated by a parts integration of the 
terms containing the singularities. After the first-order infinity has 
been treated and the x^_-, x-, and z-lengths nondLmenslonalized, equa- 
tion ( 23 ) is 


.7 
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Numerical calculations may be made by using equation (2b) for all 
values of 0 O except 9 0 = 1 where f Q becomes indeterminate and. the 

arc hyperbolic functions become infinite. The indeterminacy when 
evaluated, yields 


v 2 o 2„ 2 

(*> \ -*0 " °o ^o 

' °'®o = l 2y~ 0 
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Thus, an integration by parts of the singular terms of equation (25) 
similar to that made to obtain equation (24), using now ( f o)e o= i 

expressions for k L and kg in which 0 O Ms been set equal to 1, 

yields 





P\I ~ e 2 ~ ^ 

0 ) kg 



The integrations in equation ( 25 ) may be handled by numerical methods. 
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Sidewash due to doublets distributed In wake in region D .- In 
region D, equation (l^) takes the form 


Xd 


Urn 
y — X> 


zHB 2 d 


2 §y 


Z°Z. 

J a ^ pyu/so 2c2 - 




2 2 
^+2 


dxi 


[jx - c) 2 - - p 2 2 2 j3/2 


(26) 


Carrying out the differentiation and then the first integration for y = 0 
gives 


9 0 c 

zH(x - c) P P 

3p 2 yi 2 J Q o 2 ° 2 - P 2 yi 2 

*p 2 Jo 

^y^ 2 + z 2 )jjx - c) 2 - p 2 y^ 2 - p 2 ^ 2 1 ^ 2 


(x - c) 2 g 2 y I 2 ye o 2 c 2 - P 2 ^ 2 

(yq 2 + z 2 ) 2 (x - c) 2 - p^yj 2 - p 2 z 2 p/ 2 



By making the variable substitution 



9q 2 c 2 

P 2 


and defining 


9q 2 c 2 

(x - c) 2 - p 2 z 2 
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equation (27) reduces to 



The expression for ^ given by equation (28) may be integrated (see 
appendix B) to yield a closed-form solution for v w p. The nondimensiona] 
v 

TJ T) 

equation for — — is 
pb/2 


t 






^o^oN ~ ^o , 
e o( z o 2 + l) 


' (29) 
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Sidewash due to doublets distributed In wake of region E .- In 
region E the derivation of the sidewash is similar to that of region D 
and yields 


V W,E _ 2 (*o ~ X ) 
ib72 *e 0 G(e 0 ) 



+ 


^-l,o^o z o' 




+ 1 


Xq - 1 


+ 




(30) 


where 


%o 


2 _ (x - e) 2 - ft 2 z 2 


0o 2 - 2 


Sidewash at x-axls . - In the xy-plane (plane of the wake), only the 
doublets distributed in the wake contribute to the sidewash. Equation (20) 


for 


P,D 

pb 72 


and equations (24) and (25) for 


approaches zero, whereas equations 


P,E 
pb/2 

(29) and (30) give 


approach zero as z f 


V W,D = V W,E = 1 

pb/2 pb/2 G(0 O ) 


( 31 ) 



NACA TN 3609 


2k 


This result is identical -with that which would he obtained by use of the 
formula obtained in reference J by considering the properties of vortex 
sheets. This formula, for y = y^ = 0, is 


v 

ib72 


(x,0,0+) 



( 52 ) 


where 



( 55 ) 


Sidewash at x = <*> (Trefftz plane) .- As x approaches infinity, the 
contribution of the doublets distributed on the plan form to the sidewash 
goes to zero and the total sidewash is given by 


V W,P\ = / y \ = 1 A + 2z Q 2 

\pb/2/ 

x — G ( 6 o)\f7T7 



(3*f> 


Equation (3*0 could also have been determined more directly by using the 
formula 



lim 
y — >0 


z 

_d_ 

2itpb/2 

Sy 



(y- - yi ) 2 + z 2 


(55) 
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which, may he obtained from equation ( 2 6 ) by performing the first integra- 
tion and then setting x equal to infinity. 


Lifting-Line Sidewash 

The lifting-surface method by which the sidewash behind a rolling 
delta wing was derived in the previous section is applicable to wings of 
arbitrary plan form; however, the integrations which would be required 
before the potential or one of the perturbation velocities could be 
obtained in a calculable form are extremely difficult to evaluate. It 
is of importance, therefore, to develop some approximate expressions 
which may be easily evaluated either analytically or numerically. Refer- 
ences 2 and 8 indicate that a lifting line and a lifting line approximated 
by supersonic horseshoe vortices can be used as good approximations to 
lifting-surface solutions for most downwash problems. It would seem that 
a comparison of the sidewash behind the rolling delta wing calculated by 
the lifting-surface method with that calculated by an approximate method 
might give some indication of the usefulness of the "approximate" approach 
for sidewash problems. 

References 2 and 8 together represent a fairly thorough study of the 
lifting-line and approximate lifting-line methods, especially with regard 
to downwash calculations, and show that the bent (swept) lifting lines 
will probably give the best results for swept and triangular wings. The 
potential due to a bent lifting line may be obtained from the errata of 
reference 2 as 



r ( y l) ^-1 zh- 2 - P 2 fe 2 + z 2 ) 


2ir 


YX 


m 


m 


1*2 




_1_ 

2rt 



ar fa) 

ayp 


tan 


-1 z 


'ft 2 


3 2 (y 2 + z 2 : 


£ 

m 


YX - — - — 


2 

m 


ay-. 


(36) 


where the equation of the lifting line is 


X 1 = 


y x + k 


m 
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and the circulation r is defined hy equation (33) 
across the surface evaluated at the trailing edge. 


*1 


is zero, equation (3 6) becomes 


as the potential jump 
(See fig. 3-) When 



m m 


hg 

h l 



which might he considered as the potential in space (at a point x, y, z) 
of a finite bent vortex of constant strength. A number of finite bent 
vortices distributed along a line can be used to approximate the potential 
in space of a lifting line with any prescribed lift distribution. 

From equation (37) , the sidewash due to a bent vortex is readily 
obtained by taking the derivative with respect to y. The following 
result is obtained: 


v 


M = JL 

dy 2 it 



(58) 


When m approaches infinity, equation (38) becomes the sidewash for a 
rectangular horseshoe vortex and agrees with the equations given in 
references 9 and 10. 

Since the loading on a rolling wing is antisymmetrical, the induced 
sidewash from each panel is in the same direction and equal in the 
y = 0 plane. For this reason it is necessary to calculate the sidewash 
only from one panel and double it. 

Equation (38) can be utilized to formulate an apprQximate expression 
for the sidewash due to a series of constant-strength horseshoe vortices 
spaced along a 1 i ne so as to represent as closely as possible some pre- 
scribed span loading. This expression is 
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r ( y i+i) ' r ( y i-i) 

feo 'w 


zY. 


^(gX^ - p^ 2 - p 2 z 2 ) - zXjfc 2 - p 2 z 2 ) 


/ X i 2 - P^i 2 ' g 2 z 2 (y* - ^ - Jjfj* + z2( Xl 2 - p 2 ^ 2 . p2 z 2) 


( 39 ) 


Y-i + k 

where = y - y±, = x - Xj, and xi = — , the subscript i 

^ x m 

takes on all integral values from 0 to n. Equation ( 39 ) in nondimen- 
sional form is 


V 

pb/2 



L( y A±i). 

p(b/2) 2 


r fa-l)' 

p(b/2) 2 




fa - P 2 Y ij0 2 - P 2 z 0 2 ) - - P 2 z 0 2 )' 




fa - 


P^Y. 


1,0 


- 3 c z, 


2 , 2 


ijO 5 ! 



(J.0) 


In the application of equation (40) to the calculation of sidewash, 
some care should be exercised that the forecone from the field point 
under consideration does not intersect the lifting line at a point close 
to the corner of a horseshoe vortex. When the forecone intersects the 
lifting line near a corner located within the forecone , the expression 
under the radical in the denominator of equation (10) becomes small and 
the sidewash becomes large. (See sketch 1 below. ) 



Sketch 1 


Sketch 2 
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A zero value for the square root and an infinite value for the sidewash 
result when the forecone intersects the corner. The abruptness of the 
Infinity varies with the distance of the field point from the corner. 
Note that when x = «>, the infinity no longer exists . The preferable 
field-point locations have forecones intersecting the lifting line as 
shown in sketch 2. The closeness of the forecone to the corner when the 
corner is outside the Mach cone is obviously of no consequence because 
only line vortices within the forecone from the field point contribute 
to the sidewash at the point. 


RESULTS AND DISCUSSION 
Exact Sidewash 


The exact linearized nondimensional sidewash has been calculated 
for values of 0 O of 1.00, 0.75> O.5O, 0.k0 y and 0.30; values of x Q 

from 1.2 to 2.4; and values of z Q from 0 to 0.6 except where these 

values are ahead of region E. Variations of the sidewash param- 

v/v s 

eter — — with Xp. for 7 values of z 0 from 0 to 0.6 and for the 
pb/2V 0 0 

0 o 's given are presented in figure 4. Cross plots of figure 4 which 
show the variation of the sidewash parameter with z 0 for 7 values of 
Xq from 1.2 to 2.4 are given as figure 5. 

In order to depict the effect of Mach number and leading -edge sweep, 
variations of the sidewash parameter with z D for values of 9 0 of 1.00, 

O.75, 0.50, and 0.30 have been plotted for three longitudinal locations: 
x Q = 1.6, x 0 = 2.0, and x Q = «> (fig. 6). (it should be noted at this 

point that an increase in 9 Q may be interpreted as either an increase 

in Mach number for a fixed leading-edge slope or an Increase in the wing 
semiapex angle for a specific Mach number.) The major difference to be 
noted in the effects of changing 9 0 is that, when the longitudinal sta- 
tion is ahead of the line of intersection of the Mach cones from the 
trailing-edge tips, an increase in 9 0 causes an increase in the sidewash 

at the higher values of z Q , whereas the sidewash at a Btation remaining 
behind the intersection line during an increase in 9 0 experiences a 
decrease in sidewash at all values of z Q which are unaffected by the 

localized infinity at the intersection line. By way of illustration it 
can be seen at station Xq = 1.6 that when 9 0 is increased from 0-75 

to 1.00 the sidewash increases for values of z Q greater than 0.2; for 
values of 9 0 of 0.3 and Q.5> when Xq = 1.6 is behind the intersection 
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line, the effect of increasing 0 O is to decrease the sidewash at all 
values of z Q except at z Q = 0.6. This point is affected by the infin- 
ity at the intersection line. 


Approximate Sidewash 

From the nature of the analytical and numerical integrations required 
to obtain the exact sidewash for triangular wings, it is apparent that for 
wings with more complex potential-jump expressions the derivation of exact 
sidewash would be a difficult task. Herein lies the merit of the approxi- 
mate lifting-line method (eq. (^)) which is not encumbered by the com- 
plexity of the wing-loading expression. The approximate method is, how- 
ever, hindered to the extent that an area distribution of loading is 
assumed concentrated on one or several lines . The penalty that this 
assunption imposes on the quality of the results cannot be ascertained in 
every case. For the triangular wing treated herein, comparisons may be 
made between the results from the lifting-line and the lifting-surface 
methods, and perhaps some indication may be obtained as to the regions 
wherein the approximate method may or may not give reliable values. 

By use of equation (40), approximate lifting-line calculations of the 
sidewash have been made for 0 O = 1.00 and 0.40, values of x Q from 1.2 

to 2.h, and values of z Q from 0.1 to 0.6. Sidewash values for z Q = 0 

were obtained from equation (32). A comparison of the sidewash calcu- 
lated by the lifting-line and lifting-surface methods has been made in 
figure 7, and the agreement is shown to be good everywhere except at the 
high values of z Q close behind the trailing edge for 0 Q = 1.00. The 

agreement in this region is considerably better at 0 O = 0.40 (fig. 7(b) ) 
and indicates that as 8 0 is decreased from 1.00 the approximate calcu- 
lations will become more reliable at locations close behind the trailing 
edge. 


Seventeen bent horseshoe vortices were used to approximate the 
lifting line, with the concentration of vortices greater near the 
tip because of the rapid change in the span loading in this region. 

The shape of the bent lifting line used in the approximate calculations 
was a pair of straight lines connecting the midpoint of the root chord 
to the tips . Additional computations of the sidewash have been made 
using lifting lines composed of straight lines connecting the. tips with 
the c/k point and connecting the tips with the 3 C /^ point but the 
agreement with the exact sidewash was not so good as that evidenced in 
figure 7 • 
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It is of interest that the spanwlse center of loading of the loading 
distributed along the lifting line connecting the c/2 point with the 
tips (the lifting line yielding the best agreement with the lifting- 
surface results) was located longitudinally closer to the actual wing 
center of loading at ^Q-Jb than it was when the loading was distributed 
on the other two lifting lines. 

The wing loading in the examples just discussed was distributed on 
one lifting line. Sidewash obtained by distributing the wing loading on 
more than one lifting line would probably show better agreement with the 
exact results in region E, because some -effect of the longitudinal dis- 
tribution of loading over the wing could then be realized. 


Example of Flow-Field Effect on Vertical. Tail 

The effect of the induced sidewash velocity behind a rolling wing 
on the forces and moments contributed by a vertical tail can best be 
illustrated by analyzing a specific wing -tail configuration. The perti- 
nent geometric characteristics of the wing-tail model are (see fig. 8): 


Wing aspect ratio .... 
Tail aspect ratio .... 
Tail area 
Wing area 
Tail span 
Wing semispan 
Tail chord 
Wing chord 

Center-of -gravity location 


5-2 

1.866 

0.21 


0.7 


0.6 


0.5c 


A free-stream Mach number of 1.6 (3 = 1.25) has been chosen. The wing 
leading edge for this Mach number is sonic (0 O = l) and the vertical- 

tail leading edge is supersonic. Induced side-force and yawing moments 
for a number of ■ longitudinal positions of the vertical tail have been 
determined by numerical integrations in a manner similar to that used in 
reference 11 to obtain the contribution of horizontal tails with super- 
sonic leading edges to the lift and pitching moment. In making the 
numerical integrations, sidewash curves (figs. ij- and 5 ) were used which 
had the infinity at the tip-cone intersection line faired through. 
Isolated vertical -tail forces and moments have been computed from the 
formulas given in reference 12. 
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Figure 9 shows in stability-derivative form the variation of the 
induced, isolated, and total forces and moments with the longitud in al 
location of the vertical tail. For the example configuration chosen, 
the induced forces and moments are greater than the "isolated" forces 
and give rise to a positive Cy^ coefficient and a negative C n ^. 

Obviously, from the sidewash curves, if the vertical tail were moved 
away from the x-axis, the induced force would be reduced. The isolated 
forces and moments, on the other hand, would increase and the total 
Cy and total Cl would become negative and positive, respectively. 

P T> 


Assumptions and Limitations 

In some cases, the assumptions made in the analysis, by necessity 
or for convenience, to allow the determination of the sidewash behind 
triangular wings by linearized supersonic flow theory limit the applica- 
tion of the results. Some discussion of these assumptions and limita- 
tions may be useful. 

The validity of the assumption of a flat vortex sheet for wings 
with very low aspect ratios is questionable, but, in the absence of 
experimental and theoretical information directly concerned with the 
vortex sheet behind rolling wings, no definite statement can be made as 
to the effects that wing aspect ratio, roll velocity, and distance behind 
the trailing edge will have on the rolling -up of the vortex sheet. It 
may be possible, as suggested in reference 13, to get some indication of 
these effects from the data published in references l4 and 15 concerning 
the rolling -up of the vortex sheet behind wings at an angle of attack. 

The vortex sheet has been assumed not to rotate. The angle (in 
degrees) through which the vortex sheet would rotate in moving from the 
wing trailing edge to a point d distance downstream is given by 

Tv _ 180 

* V 


= 180 pb ( d \ 

* 2V \b/2 / 
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Substitution into this formula of values of pb/2V usually encountered 
and of distance d up to two semispans will generally yield rotation 
angles small enough to be neglected. 

The sidewash expressions derived in the analysis for rolling wings 
at zero angle of attack are also applicable for finite angles of attack 
because angle -of -attack loadings are symmetrical and do not contribute 
to the sidewash in the xz-plane , The displacement of the assumed flat 
vortex sheet from its zero angle -of -attack position, however, must be 
accounted for; that is, the sidewash given for a point x, z=0 + for the 
zero angle -of -attack case represents the sidewash at the vortex sheet 
when the wing is at an angle of attack. At a distance d behind the 
trailing edge, the displacement of the vortex sheet below the trailing 
edge may be found (see following sketch) from 


Id 

h = tan e dx 

UT.E. 


^Wing 



Values of tan e for a triangular wing are given in reference 2. 


CONCLUDING REMARKS 


The variation of sidewash with longitudinal distance in the verti- 
cal plane of symmetry behind rolling triangular wings traveling at super- 
sonic speeds baa been derived by linearized lifting-surface and lifting- 
line methods. The range of supersonic Mach numbers for which the 
liftin g -surface results are valid is limited by the condition that the 
wing leading edges must be subsonic. The variations of lifting-surface 
sidewash are presented in graphical form for a number of values of 0 O , 
a Mach n umb er — leading-edge-sweep parameter. Sidewash calculated from 
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the lifting-line formula has been compared with the lifting-surface side- 
wash for values of 0 O of O.hO and 1.00. This comparison shows very good 

agreement of the lifting-line results with lifting-surface results except 
at the higher vertical distances close behind the trailing edge for 
0 O = 1.00. The curves for 0 O = 0.40 reveal that, as 0 O is decreased 

from 1.00, the agreement close behind the trailing edge improves. 

An illustrative calculation of the sidewash-induced force on a half- 
delta vertical tail operating behind a triangular wing indicates that the 
induced force acts in opposition and is comparable in magnitude to the 
damping force created on the isolated rolling tail. In order to determine 
the total force which would act on the vertical tail of a - steady-rolling 
aircraft in flight, an additional force induced by the aileron sidewash 
should be calculated. This force may be of the same order of magnitude 
as, and opposed to, the force induced on the tail by the W^Lng sidewash. 

No attempt has been made in the present paper to evaluate aileron sidewash. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., October 21, 1955 • 



3k 


NACA TN 3609 


APPENDIX A 

TREATMENT OF FIRST-ORDER INFINITIES IN Vp^-g 
AS GIVEN BY EQUATION (20) 


In order to isolate the "infinite" term, equation (23) may he 
written as 



cbci 


© 


&2r 


E r 


Jo 6 o x l 1 - 




dxp 

®o x l 



+ 





dxp Ep 

e o x l 1 - k x 2 


(Al) 


Terms (l) and (3) are integrahle hy numerical methods, whereas terms 
(2) and (h) contain a first-order infinity at the limit xp = f 

(when Xp = f , kp = kp = 1^ . 

Consider term ( 2 ) in equation (AL) 



dxp kgEg 

e o x l 1 - kg 2 
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If k 2 is replaced by its equivalent (see symbols), term (2) becomes 


0 (x - Xp) 2 - 0 2 Z 2 - Qq^Xj 2 


(x - xj^ - p 


2 - S 2 z 2 


Assume ejqpression (A2) to- be of the form / u dv, where 

do 


u = E 2 /(x - xj 2 - p 2 z 2 


(x - xj 2 - p 2 z 2 - Q 0 2 x 2 


Then, 


coth' 


_1 - (l - Bp 2 )*! + x 


jq q 2 x 2 + p 2 z 2 (i - e 0 2 ) /e/x 2 + p 2 z 2 (i - e 0 2 ) 


e o( E 2 “ K 2) K 2 k 2( x - x^) 


e o x l 
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and a parts integration of expression (A2) gives 


. : x 0 2 - p2 . z ?_ coth -i . 

- i 1 ~ 0 o 2 ) x i + x 

ye 0 2 x 2 + P 2 z 2 (1 - 0 O 2 ) ^ 

/e 0 2 x 2 + p 2 z 2 (1 - e Q 2 ) 


0 


coth - - 

L _ 

-I 

CM 

0 

<33 

1 

<H^ 

lx-j_ + x 



J 

'e/x 2 + p2 z £ 

! (^ - »o 8 j 

0 o(E2 ~ 


£ 

2x2 

+ P 2 Z 2 ^ 

1 

CD 

O 

1 

kg 


~ x l) 
e o x l 


<2^ (A3) 


Substituting for its equivalent in term (j+) of equation (Al) results 

in 



9 o x l E i toj 

- (x - x-^) 2 + 0 2 z 2 


(m 


With 


u = QqX^ 


and 


dv 


-dx-^ 

Qq 2 *! 2 _ ( x ~ x l) 2 + P 2 z 2 
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a parts integration of expression (A4-) yields 



The nonintegral term in expression (A5) may "be written in slightly 
different form as 



Inasmuch as when Xp = f, kp = kg = 1^ it is clear that the evaluation 

of the integrated term in expression (A5) (equivalent to expression (a 6)) 
at the limit xp = f will cancel the integrated term of expression (A3) 

evaluated at this same limit. The complete expression for Vp^-g is now 
seen to he 
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This equation in nondimensional form is presented in the analysis of the 
report as equation (24). 
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APPENDIX B 


EVALUATION OF v w ^ D 


The contribution of the doublets distributed in the wake to the 
sidewash in region D is given by equation ( 28 ) as 



which by use of partial fractions can also be written 
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For ease in writing, -2— has "been replaced "by a in equation (B2). 

zp 

Consider first the integral 




ds 


(l + a 2 s 2 )/l _ ^ 


2 s 2 


(B3) 


If the variable transformation s = sin 0 is made, expression (B3) 
becomes 



cos 2 0 d0 


fl + a 2 sin 2 0^ Jl-k^^ 2 sin 2 0 


(B4) 


The evaluation of expression (B4-) is given by formula (9), table 6l of 
reference l6 as 



A 


s 2 ds 


(l + a2s2 )i 1 ~ k 2,o 2s2 


a/ 8,2 + k 2,o 2 


(B5) 


where 



|+ (k ^ 0 - Eg^jF^ot-ia,!^^*) - ^^(cot^a,^^*) 


(b6) 



6T 

NACA TN 3609 41 


Expression (B4) could also have been integrated without recourse to 
tables with the aid of reference 17 (pp. 134 to 136) . The integration 
of expression (B3) 



may be performed by using the relationship (see p. 13 of ref. 17 and 
p . 79 of ref . 18 ) : 



From equations (B5) and (b 6), equation (B7) becomes 
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Carrying out the differentiation in equation (B8) results in 



\Jl - s 2 ds 

(l + - k2 j0 2 a 2 


«V a + ^.o 


1 

2 


ava4 + 1 



+ 


l /(* 2 + i) (a 2 + k 2 j0 2 ) 


G ( k 2,o> a ) + 


1 

2 



8,2 + ^.Q 2 ^,0 ~ E 2,o 

(a 2 + l) 3 ^(a 2 + l) (a 2 + k^ 0 2 ) 


(B9) 


The remaining integration needed to evaluate 


W,D 


(eq. (B2) ) is 



JT 


s 2 ds 




k 2 s 2 

K 2,o 8 


This integration may he reduced to standard elliptic forms hy the 
transformation s = sn u and integrated to give 


f 1 _ £ZZ- §i % ? E ? ,° (bio) 

0 ( l-^,o Z ° 2 Y /2 


The side-wash ^ is now completely defined hy equations (B2), (B5), 
(B9), and (BIO) as 
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and 




Hpe 0 c 

2 + a 2 £ 


ay? + 1 


a(x - c)}^ 

f (k2 ’°’ a) ' e 0 c( a 2 T^ 2 ’' 


Replacing a by its nonfli mensional equivalent 


a 


9 o c _ _1_ 
pz Z Q 


and H by its equivalent and then nondimensionalizing x by 
v w D by pb/2 gives 


V W,D _ 2 

2 *o + 1 

L _i) 

K2,o^2, 

o( x o - ^) z o 

pb/2 itG(0 o ) 

N 

O 

ro 

+ 

H 


9 o 

(*° 2 + x ) 


Equation (B13) is identical with equation (29) of 
the exception that, in the analysis, the function 


the analysis 


(B 12 ) 


c and 


(B13) 


, with 
has been 


written out. 
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(e) 0 O = 0.30. 

.Figure 4 .- Concluded. 
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(a) e 0 = o.4o. 


Figure 5*~ Continued. 
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Figure 5*- Concluded. 
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Figure 6.- Concluded. 
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Figure 8.- Wing-tail model used to ill u strate effect of wing sidevash on 

a vertical tail. 
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